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ABSTRACT 

We present a small computer algebra program for use in Maxwell's theory. 
CO ' The Maxwell equations and the energy-momentum current of the electromag- 

netic field are formulated in the language of exterior differential forms. The cor- 
. responding program can be applied (in the presence of a gravitational field) in 

curved Riemannian spacetime as well as in flat Minkowski spacetime in inertial 



in 

On 



or non-inertial frames. Our program is written for the computer algebra system 
REDUCE with the help of the EXCALC package for exterior differential forms. 

The two major advantages of this modern approch to electrodynamics — the 
natural formulation, free of both metric and coordinates, and the straightforward 
bJ[). programming of problems — are illustrated by examining a number of examples 

ranging from the Coulomb field of a static point charge to the Kerr-Newman 
^ ■ solution of the Einstein-Maxwell equations. 

5h 

1. Introduction 

In the lecture of Hartley!, the effective use of exterior differential forms has been 
demonstrated convincingly. We were led to our Maxwell program because we were 
searching for electrically charged solutions in the framework of general relativity — see 
the talk of Brans! in this context — and of alternative gravitational theories!, and only 
streamlining our program using exterior differential forms made it a feasible task. We 
believe that our program, besides being helpful in physics, could also have applications 
in electrical engineering (see the textbook of Meetz and Engl!, written for engineering 
students, inter alia, where the calculus of exterior forms is applied throughout). 

*Supported by the Graduiertenkolleg 'Scientific Computing', Koln-St.Augustin. 

t Computer Algebra in Science and Engineering. Proc. of the ZiF- Workshop, Bielefeld, Aug. 1994. J. 
Fleischer et al. (eds.) World Scientific, Singapore, to be published (1995). 
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Fig. 1. Graphic representation of geometric objects used in Maxwell's theory (adapted from 
Burke!'! and Konzeni). 



Our axiomatic formulation of Maxwell's theory is based on three principles stand- 
ing on firm experimental bases: (i) the conservation of electric charge, (ii) the notion 
of the Lorentz force, and (iii) the conservation of electromagnetic flux. 

2. Maxwell's Theory in Brief 

2.1. Electric Charge Conservation 

The inhomogeneous set of Maxwell's equations, in terms of exterior differential 
forms, may be written as 



The sources are the electric current density 2-form j and the electric charge density 
3-form p. The magnetic excitation 1-form is denoted by 7i and the electric excitation 
2-form by T>. The underline d is the 3-dimensional exterior derivative, and the dot 
denotes the time derivative. 

These equations, reformulated in 4 dimensions, read 



provided we introduce the time-space decompositions of the electromagnetic excitation 
2-form (t = time coordinate) 




(1) 



dH = J, 



(2) 



H 



HAdt + V 



(3) 



and the electric current 3-form 



J 



j Adt + p. 



(4) 
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If the spacetime under consideration is contractable, then it can be shown that (|]) can 
be derived from the integral version 



J= / dJ = (5) 

dV 4 Jv A 



of the conservation law for electric chargemld. 

2.2. Lorentz Force 

Postulating the (vectorial) Lorentz-force equation as^] 

f = p A E + j A B , (6) 

we are led in a most natural way to the electrical field strength 1-form E and the 
magnetic field strength 2-form B, with E = e-iE and B = e_i£>, where e = {ex, e 2 , e 3 } 
represents an arbitrary 3-dimensional vector basis (frame), a so-called triad, and _i 
denotes the interior product. In analogy to (|3|), the field strength can be collected in 
the 4-dimensional electromagnetic field strength 2-form 

F = —dt A E + B . (7) 

One possible graphic representation of the sources and fields of Maxwell's theory is 
displayed in Fig.p]. We find this way of representing exterior differential forms quite 
intuitive and, for more details, refer to the literature cited in the caption of Fig.|l]. 

2.3. Electromagnetic Flux Conservation 

The homogeneous set of Maxwell's equations 



dE + B = , 
dB = 0, 



(8) 



can be rewritten by means of ([?D as 

dF = , (9) 
or, in integral form, as an electromagnetic flux conservation law 

F= [ dF = 0. (10) 

£>V 3 Jv 3 



*In 4 dimensions, instead of (||), we have f a = (e a -lF) A J, with a 4-dimensional frame e a 
{e ,ei,e 2 ,e 3 }. 
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2.4- Maxwellian Structure and Spacetime 

So far we have the two Maxwell equations (||) and (fj|) and the expression for 
the Lorentz force (||) interrelating the current and field strength by means of the force 
concept of mechanics. These equations can be formulated on any 4-dimensional differ- 
ential manifold, as long as it allows a (l+3)-decomposition. It is worth noting that 
these Maxwellian structures emerge prior to the introduction of a metric or a connec- 
tion. The conservation laws (||) and (|1(]) only require the possiblity to circumscribe a 
boundary around an arbitarary 4- or 3-dimensional volume element, and to count the 
charge or flux units contained therein. 

In particular this means that the equations (|2]) , ([|) , and ((^) are valid in this form 
in a curved and contorted spacetime, i.e. with non-vanishing curvature and torsion, or 
in a flat Minkowskian spacetime in inertial and non-inertial frames - they always keep 
this form. 

However, in (|2|) and (|9|) we have only 6 independent evolution equations, namely 
(HDi and for specifying the 2x6 independent fields {H,F}. Something is still 
missing, namely the constitutive laws, such as the one relating H to F, and Ohm's law. 

In a vacuum, the field strength and excitation are related by 

1 *ti .„ ,. action V 
H=-*F, with dunu = — - = — = £!, (11) 

fj, charge 2 A 

where * is the Hodge star operator, which does depend on the metric — or rather on 
its determinant. 

3. Program 

In an interactive Reduce 3.5 sessionEiHll, the Excalc packaged can be loaded by 
the commancfj] 

load_package excalc$ 

First Excalc has to be acquainted with the rank of the differential forms to be 
used in the code. This is achieved by the declaration 

pform potl=l, {f arad2,excit2}-=2, {maxhom3,maxinh3}=3$ 

where potl and f arad2 are the identifiers of the electromagnetic potential and field 
strength, excit2 represents the electromagnetic excitation, and maxhom3 and maxinh3 
denote the left hand sides of the homogeneous and the inhomogeneous Maxwell equa- 
tions. The numbers after the equal sign specify the rank of the forms. 

The mathematical symbols used in exterior calculus, translate into Excalc as 
listed in table [l]. 

has written the Mathematica package CARTAN which should be also useful for program- 
ming Maxwell's theory. 
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Math. Symbol 


Excalc Symbol 


Operator 


Operator Type 


A 

A 




exterior product 


"nary" infix 


d 




p^rtprior npriA/nfivp 

CA Id IvI VJ-VT/X 1 V CL u 1 V 




d 


Q 


partial derivative 


"nary" prefix 


_i 


J 


interior product 


binary infix 


£ 


1- 


Lie derivative 


binary infix 


* 


# 


Hodge star operator 


unary prefix 



Table 1. Translation of mathematical symbols into Excalc. 

Now we are ready to write a small prototype program. Suppose we consider 
Minkowski spacetime and use polar coordinates {x°, x , x 2 , x 3 } = {t, r, 9, </?}. Then we 
have to specify a local basis of 1-forms, a coframe, by means of the Excalc declaration 

coframe o(t) = d t, 

o(r) = d r, 

o(theta) = r * d theta, 

o(phi) = r * sin(theta) * d phi 
with signature (1,-1,-1,-1); 
frame e; 

The command frame e; provides us with the frame e a dual to the coframe that we 
have chosen to call o a . 

3.1. Static Point Charge 

As test scenario, let us take the Coulomb field of a static point charge q. The 
corresponding electromagnetic potential 1-form A and the field strength 2-form F = 
dA in Excalc read 

potl := -(q/r) * d t; 
farad2 := d potl; 

Now we translate into Excalc the homogeneous Maxwell equation (|9|) 

maxhom3 := d farad2; 
the constitutive law for vacuum ([TT| ) 

excit2 := (1/mu) * # farad2; 
and the inhomogeneous Maxwell equation (^) 

maxinh3 := d excit2; 
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The result of these commands, typed in during an interactive Reduce session, is that 
both left hand sides of Maxwell's equations are zero; in other words, Maxwell's equa- 
tions are fulfilled, and the ansatz for A is justified. 

The canonical energy-momentum 3-form of Maxwell's field is given by^] 

E a = e a jI + (e ft jF)Ai?, (12) 

where L is the Lagrange 4-form of Maxwell's field: 

L = ~FAH. (13) 

We declare the corresponding forms 

pform lmax4=4, maxenergy3(a)=3$ 
and define the Lagrangian 

lmax4 := -(1/2) * farad2 ~ excit2; 



Now the energy-momentum 3-form (|T2|) can be computed straightforwardly (note that 
a lower index is always marked with a minus sign in Excalc): 

maxenergy3(-a) := e(-a) _| lmax4 + (e(-a) _| farad2) ~ excit2; 

We find 

E t = -^o r Ao e Aof (14) 

2/ir 4 

and (up to minus signs depending on the ordering of the basis 1 -forms o a ) equivalent 
expressions for the three other components. As expected, energy-momentum turns out 
to be proportional to q 2 /r 4 . 

3.2. Small Ring Current 

An equally simple task is the potential of the magnetic dipole moment A4 of 
a small ring current. From the literature (see Jacksonll3, p. 182) we take the vector 
potential in the Minkowski spacetime as A(r) — M. x r/r 3 . We evaluate the vector 
product (see Fig. |2|) and find A = Ai/r 2 sin 9 e v , where e v is the unit vector in <p- 
direction. Then the four-potential 1-form reads 



A = ^ sin 9 o v = — sin 2 9 dp. ( I V) 



* The conventions we use for energy-momentum are as follows: define the energy-momentum tensor T(j a 
in terms the energy-momentum 3-form E Q by S Q = Tp a rf with t/ 3 := *o . If this Tp a is multiplied 
by —1, we have the conventions of Ref.16. This is the reason for the minus sign on the right hand side 
of (0) below. 
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Fig. 2. The electromagnetic field of a small ring current. 



Let us denote the angular momentum of the electrically charged circling particles 
by L and the angular momentum per unit mass by a := L/m. Then the gyromagnetic 
parallelism yields (Jackson^ p. 183) M. = gqa/2 with q = charge or, substituted in 

©, 

A = sin 2 6 dep. (16) 
2r 

By g we denote the ^-factor of the current distribution, which equals 1 for orbital 
angular momentum (as for a ring current) and 2 for spin angular momentum. 
In Excalc we have 

potl := g*q*a/(2*r) * sin(theta) **2 * d phi; 
and, for field strength and excitation, 



farad2 
maxhom3 
excit2 
maxinh3 



= d potl; 

= d farad2; 

= (1/mu) * # farad2; 

= d excit2; 



The Maxwell equations turn out to be fulfilled. We will apply the potential (jig) in 
Sec.O in a less trivial context. 

4. Examples 

4-1. Electromagnetism in a Rotating Frame 

In the flat spacetime of special relativity, inertial frames of reference are realized 
by (pseudo-)orthonormal tetrads which are all arranged in a parallel fashion — here 
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parallel is meant in the ordinary (pseudo) Euclidean way. A corresponding coframe, in 
polar coordinates, was specified in Sec.|3[ If we want to describe non-inertial frames of 
reference, then we can turn to an accelerated and rotating frame, as used in Ref.18, for 
example, (see also Ref.19), and the parallel arrangement of the frames will be disturbed 
thereby. 

There is a plethora of different books and papers on electrodynamics referred to 
non-inertial frames, see Refs.20,21 and literature given there. Formulated in terms of 
exterior differential forms, such formalisms simplify appreciably. We will demonstrate 
this claim by investigating the question of how a rotating observer will see the field of 



the magnetic dipole as specified in Sec.pT2 



Let us assume that the axis of rotation coincides with the z-axis such that lo is 
the only nonvanishing component of the local angular velocity. We use the coframe^] 
as given in ThirringEHl (last equation on p. 556) in order to express this configuration in 

Excalcfl: 

coframe 

o(t) = (d t - omega*rho**2* d phi)/sqrt(l - (omega*rho) **2) , 
o(rho) = d rho, 

o(phi) = (d phi - omega* d t) *rho/sqrt (1 - (omega*rho) **2) , 
o(z) = d z 
with signature 1,-1,-1,-1; 
frame e; 

The dipole potential (|H]) is expressed in terms of polar coordinates. To put it into our 
rotating cylindrical coordinates, we would need a pullback operator for forms which 
will, however, only be available in the next Excalc release. We shall overcome this 
deficiency by a small trick. Define 

pform {r ,theta}=0$ 



r := sqrt(rho**2 + z**2) ; 

theta := acos(z/sqrt(rho**2 + z**2)); 

If we start with polar coordinates {t, r, 8, if}, Excalc will then substitute cylindrical 
coordinates {t, p, <p, z} instead. Using this transformation, the dipole potential ( |16|) 
can be viewed from the rotating system of reference: 

factor o(0) ,o(l) ,o(2) ,o(3) , ~$ off exp$ on gcd$ 
pform pot 1=1, {f arad2,excit2}=2, {maxhom3,maxinh3}=3$ 

potl := g*q*a/(2*r) * sin(theta) **2 * d phi; 
farad2 := d potl; 

§By using Excalc, one can prove easily that the Riemannian curvature belonging to our new coframe 
vanishes identically, that is, we are still in a flat Minkowski space. 

^Before making a new coframe statement, we recommend starting a new Reduce session. 
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By making use of the factor command we caused Reduce to "factor out" the terms 
given. Additionally, we turned off and on the two Reduce switches exp (expand) and gcd 
(greatest common divisor), respectively, in order to simplify the resulting expressions. 

Let us have a closer look at the physical situation as seen by the rotating ob- 
server. For this purpose, we compute the components of the electric and the magnetic 
field, respectively (see Ref.19 for the (3 + l)-decomposition of the electromagnetic field 
strength 2-form): 

indexrange {i, j ,k}={rho,phi,z}-$ 
pform {ee(i) ,bb(i, j)}=0$ 
index_symmetries bb(i,j) : antisymmetric$ 

ee(-i) := e(-i) _| (e(-t) _| f arad2) ; 
bb(-i.-j) := e(-i) _| (e(-j) _| f arad2) ; 

The indexrange declaration allows to have certain indices not ranging over all space- 
time dimensions but over a subset that can be chosen arbitrarily. Here we declared the 
indices and k to cover only the spatial section of spacetime. 

The rotating observer still perceives, up to a Lorentz factor 7=1/ ^Jl — (cup) 2 , 
the magnetic dipole field 

7 agq ?>pz 3pz 
= — . , ■ = -7 M- 



2 (p 2 + z 2 fl 2 " (p 2 + z 2 fl 2 



B zp = 



7 agq (p 2 - 2z 2 ) (p 2 -2z 2 ) 

Bp<t> = — } _ 9 , ..9NS/9 = ^ • A/ t- 



2 (p 2 + z 2 fl 2 ' (p 2 + z 2 fl 2 

(we substituted again the magnetic dipole moment M. = agq/ 2), but in addition he 
observes an electric field 

7 aggujp (p 2 - 2z 2 ) _ up (p 2 - 2z 2 ) 
p 2 (p 2 + z 2 f/ 2 7 (p 2 + z 2 f/ 2 

E<p = 

^ 7 aqq 3up 2 z , . 3up 2 z 

E z = ' , , , ox,/o =1 M- 



2 (p 2 + z 2 fl 2 ' (p 2 + z 2 fl 2 

whose components correspond to those of the magnetic field rotated by rr/2 around an 
axis perpendicular to the z-axis and multiplied by uo p. 

For completeness, we verify that the inhomogeneous Maxwell equation is still 
satisfied: 

excit2 := (1/mu) * # farad2; 
maxinh3 := d excit2; 
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4-2. Deriving the Reissner-Nordstrom from the Schwarzschild Solution 



To give a less trivial example, we will derive the spacetime of a spherically sym- 
metric mass distribution carrying a net electric charge q. In general relativity, the 
corresponding solution of the Einstein equation is called the Reissner-Nordstrom solu- 
tion. 

We start from the uncharged spherically symmetric solution (the Schwarzschild 
solution). Its spacetime geometry is determined by the coframe that can immediately 
be expressed in Excalc: 

pform psi=0$ f domain psi=psi(r)$ 
coframe o(t) = psi * d t, 

o(r) = (1/psi) * d r, 

o(theta) = r * d theta, 

o(phi) = r * sin (theta) * d phi 
with signature (1,-1,-1,-1)$ 
frame e$ 

Note that functions have to be declared as 0-forms also. The fdomain declaration 
specifies the independent variable(s) of a given form. The Schwarzschild function ip has 



the form ip = y 1 — 2m/r, where m is the mass of the object. It can be shown, see 
below, that this orthonormal coframe satisfies the vacuum Einstein equation indeed. 

Now we want to generalize the Schwarzschild solution to a corresponding elec- 
trically charged solution. We proceed in exactly the same manner as in Sec. [3. If the 
Schwarzschild coframe, with arbitrary ip, is already prescribed (and has been checked), 
we type in (or read from a corresponding file) the commands for the point charge q: 

pform potl=l, {f arad2,excit2}=2, {maxhom3,maxinh3}=3$ 



potl := -(q/r) * d t; 

f arad2 := d potl ; 

maxhom3 := d farad2; 

excit2 := (1/mu) * # farad2; 

maxinh3 := d excit2; 

We find that Maxwell's equations are still fulfilled, even though we did not even specify 
the function ip\ Note that this only turns out to be true since we specified the potential 
potl in terms of dt and not o(t). However, the electromagnetic energy-momentum 
3-form E Q given in ([12]) acts as a source in the Einstein equation: 

G a = -^ a . (17) 

The Einstein 3-form G a (called einstein3 in the code) is defined according to (see 
Ref.23, e.g.): 

G a := (1/2) A *{o a A 0^ A o 7 ) . (18) 
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Basically we had put a point charge on a prescribed Schwarzschild geometry — and, up 
to now, we forgot the gravitating character of the electromagnetic field. The obvious 
consequence is that the corresponding Einstein-Maxwell equation is not satisfied. In 
order to fix this, we have to "relax" the Schwarzschild spacetime; this is done by 
reformulating the coframe. Since the charged solution sought should remain spherically 
symmetric, we try to change the function if) alone. The most straightforward ansatz 
that comes to mind is 

pform unknown=0$ f domain unknown=unknown(r) $ 
psi := sqrt(l - 2 * m/r + unknown); 

The new function unknown must be determined from the Einstein-Maxwell equation 



We first have to calculate the Riemann or Levi-Civita connection 1-form T a ° 
and the curvature 2-form R 13 " 7 (with identifiers chrisl and curv2). The corresponding 



einstein3(-a) := (1/2) * curv2(b,c) ~ #(o(-a)~o(-b)~o(-c)) ; 

The Excalc command riemannconx <identif ier>$ initiates the calculation of the 
connection 1-form and stores it in identifier. If unknown vanishes, we have G a = 0, 
which proves our earlier claim that the original coframe satisfies the Einstein vacuum 
equation. 

Subsequently, we compute the Maxwell energy-momentum, see flT2|): 
pform lmax4=4, maxenergy3(a)=3$ 

lmax4 := -(1/2) * farad2 ~ excit2; 

maxenergy3(-a) := e(-a) _| lmax4 + (e(-a) _| f arad2) ~excit2; 

If the vacuum Einstein-Maxwell equation is to be satisfied, the sum of the Einstein 
3-form and (a constant times) the Maxwell energy-momentum 3-form must vanish: 

pform ode3(a)=3$ 

ode3(a) := einstein3(a) + kappa * maxenergy3(a) 

II Note that the old anti symmetric <identif ier> ; declaration is unsupported already in Reduce 3.5. 
It has to be replaced by index_symmetries <identif ier>(a,b) : antisymmetric ; 
**The transposition according to chrisl (a,b) : = christl (b , a)$ is necessary to convert the Excalc 
conventions for the connection 1-form into our own convention, see Ref.24. 



Excalc code read: 




pform chrisl(a,b)=l, curv2(a,b)=2, einstein3(a)=3$ 
index_symmetries chrisl(a,b), curv2(a,b) : antisymmetric$ 
riemannconx christl$ chrisl(a,b) := christl (b,a)$ 



curv2(a,b) 



:= d chrisl(a,b) + chrisl(-c,b) " chrisl(a,c); 
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A typical component of the output of ode3 reads 



ode3 tp 




(19) 



Requiring the coefficients to vanish yields a second order differential equation for the 
function unknown. It can be solved by means of the ansatz unknown := a * r**n; A 
short Reduce calculation yields^ 

unknown := kappa/ (2*mu) * (q/r)**2; 

This reproduces, when ip is evaluated, the well known Reissner-Nordstrom function: 

psi := psi; 



It is remarkable that, although we are now working on a curved spacetime manifold, 
not a single command of the Maxwell code had to be changed. This should clearly 
demonstrate the simplifying value of the formulation of Maxwell's theory presented 
here. 

4-3. Deriving the Kerr-Newman from the Kerr Solution 

In analogy to Sec. |4.2j , we will determine the solution to the Einstein-Maxwell 
equation of a rotating axial symmetric mass distribution with net charge q. It is called 
the Kerr-Newman solution. 

Again we take as starting point the spacetime of the corresponding uncharged so- 
lution, the Kerr solution. Using Boyer-Lindquist coordinates {t, p, 8, <p}, the orthonor- 
mal coframe reads: 

pform {rr ,delsqrt}=0$ 

f domain rr=rr (rho , theta) , delsqrt=delsqrt (rho) $ 
coframe 

o(t) = (delsqrt/rr) *(d t - (a*sin(theta)**2) * d phi), 

o(rho) = (rr/delsqrt) * d rho, 

o (theta) = rr * d theta, 

o(phi) = 1/rr * sin(theta) *(a*d t - (rho**2 + a**2) * d phi) 
with signature 1,-1,-1,-1$ 
frame e$ 

ttif ode3 was of a more complicated nature, we could feed it into Maple or Mathematica, which have 
at the moment a more complete ODE solver than Reduce. 




(20) 
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The functions delsqrt and rr depend on the parameters m (mass) and a (angular 
momentum per unit mass). They have the form delsqrt = y/ p 2 + a 2 — 2 m p and 



rr = yp 2 + (a cos6) 2 . We can prove (in the same way as in Sec. f4.2| ) that this coframe 
satisfies the vacuum Einstein equation. 

To determine the Kerr-Newman solution, we first have to specify the electromag- 
netic potential. To start with, we take the potential of a rotating charge distribution 
which we approximate to lowest order of 1/p by a static point charge and a magnetic 



dipole (see Secs.fO and K 



pform potl=l, {f arad2,excit2}-=2, {maxhom3,maxinh3}=3$ 
potl := -(q/rho) * (d t - g/2*a*sin(theta)**2 * d phi); 

Here the ^-factor is to be determined. 

In Minkowski spacetime, the homogeneous and inhomogeneous Maxwell equations 
are satisfied, as we have seen in Secs.[0] and |3]^ (a superposition of solutions is again a 
solution). However, now we have to compute Maxwell's equations on the curved Kerr 
spacetime: 

let -[(sin ~x)**2 => 1 - (cos x)**2}$ factor ~$ 



farad2 
maxhom3 
excit2 
maxinh3 



= d potl ; 

= d farad2; 

= (1/mu) * # farad2; 

= d excit2; 



We find that the Hodge star operator spoils (by means of detg) the inhomogeneous 
Maxwell equation. Since we expect that this failure is connected to the 1/p behavior 
of the potential, we change the potential by the following ansatz 

pform f=0$ f domain f = f (rho,theta)$ 

potl := -q/(rho * f) * (d t - g/2*a*sin(theta) **2 * d phi); 

Before we repeat the calculation of Maxwell's equation, we again make use of the 
Reduce switches exp and gcd: 

off exp$ on gcd$ 



farad2 
maxhom3 
excit2 
maxinh3 



= d potl ; 

= d farad2; 

= (1/mu) * # farad2; 

= d excit2; 



The expression for maxinh3 looks rather complicated — if we require it to vanish, we 
get a system of second order PDEs in two variables. We try solve it by inspection; in a 
first step first let us specify the ^-factor. There are two values for g that seem plausible: 



g = 1 corresponds to a small ring current as considered in Sec. |3.2j , and g = 2 is the 



g-factor of an electron. As maxinh3 contains (g — 2) several times as a factor, the choice 
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g := 2$ 

seems promising. Indeed, if we compute the inhomogeneous Maxwell equation again 

maxinh3 := d excit2; 

the coefficient of o* A o p A o e turns out to be quite compact. If it is to vanish, we have 
to solve 

2<?a f a df . df\ , . 

cos#— — psmv— = 0. (21) 



Solving this equation by hand, fixes the function / up to a constant (which we can 
absorb into the charge q). We substitute the solution and check the inhomogeneous 
field equation: 

f := 1 + (a/rho * cos (theta) ) **2$ 
maxinh3 ; 

Note that the inhomogeneous Maxwell equation is satisfied irrespective of the explicit 
form of the functions rr and delsqrt. 

With this result, the final form of the potential can be checked: 

potl := potl; 

Q O XT , , . 

potl := — o*. (22) 

P (a 2 cos 2 9 + p 2 ) delsqrt 1 ' 

What remains to be found is the solution of the Einstein-Maxwell equation. This 
is done in much the same manner as in Sec. f4.2| by modifying the function delsqrt. In 
addition, we now specify rr without modification: 

pform ul=0$ f domain ul=ul(r, theta) $ 

delsqrt := sqrt(rho**2 + a**2 - 2*m*rho + ul)$ 

rr := sqrt(rho**2 + (a * cos (theta) )**2)$ 

The field equation must determine the new function ul. Thus we compute the left hand 
side (the Einstein 3-form) 

riemannconx chrisl$ chrisl(a,b) := chrisl(b,a)$ 
pform curv2(a,b)=2, einstein3(a)=3$ 
index_symmetries curv2(a,b) : antisymmetric$ 

curv2(a,b) := d chrisl(a,b) + chrisl(-c,b) " chrisl(a,c); 
einstein3(-a) := (1/2) * curv2(b,c) ~ #(o(-a)"o(-b)"o(-c)) ; 

and the right hand side (the electromagnetic energy-momentum 3-form) 

pform lmax4=4, maxenergy3(a)=3$ 
lmax4 := -(1/2) * farad2 " excit2; 

maxenergy3(-a) := e(-a) _| lmax4 + (e(-a) _| f arad2) "excit2; 
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of the field equation. These yield a differential equation for the function ul: 



pform ode3(a)=3$ 

ode3(a) := einstein3(a) + kappa * maxenergy3(a) ; 

This system can be solved most easily by inspection of the (^-component of ode3: 

, acos^ dui a ko 2 — 2iiui , . R , , 

ode^ = -± (f A o A o v + — J—i o A o p A o e 23 

2delsqrt 39 2rr 4 v ; 

The first term demands that u% must not dependent on the coordinate 9. The second 
term vanishes if we substitute 

ul := kappa / (2*mu) * q**2$ 

Finally, we display the explicit form that delsqrt evaluates to: 

delsqrt := delsqrt; 



delsqrt := J p 2 - 2mp + — + a 2 (24) 
V 2/i 

Again, although we modified the spacetime geometry via the function delsqrt, we do 
not worry about the inhomogeneous Maxwell equation because it has been found to be 
satisfied irrespective of the explicit form of the functions contained in the coframe. 

This completes our deduction of the Kerr-Newman solution from the Kerr solution 
and the Einstein-Maxwell equations. This derivation is certainly not very original. 
However, the charged version of the Kerr solution with dynamic torsion! was found by 
this type of technique. 
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Appendix: Maxwell program 

% file mustermax.exi, 1995-02-15 



load_package excalc$ 



°/„ the Maxwell program is preceded by coframe and frame commands: 
% specify the coframe by filling in the exact form of the 
% appropriate coframe with your favorite editor, specify possibly 
% unknown functions (zero-forms) psi etc. and their domains 



pform psi=0$ 

f domain psi=psi(r)$ 
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coframe o(0) 



d t, 
d r, 



o(l) 
o(2) 
o(3) 



= r * d theta, 

= r * sin (theta) * d phi 



with signature (-1,1,1,1)$ 
frame e$ 

7„ here begins the Maxwell program proper: unknown functions aaO, 
% aal, etc. are prescribed for the components of the potential 
% 1-form; field strength farad2, excitation excit2, and the left 
% hand sides of the Maxwell equations are defined 

pform {aaO , aal , aa2 , aa3}=0 , pot 1=1, {f arad2,excit2}=2, 
{maxhom3 , maxinh3}=3$ 

potl := aaO*o(0) + aal*o(l) + aa2*o(2) + aa3*o(3)$ 

farad2 := d potl; 

maxhom3 := d farad2; 

excit2 := (1/mu) * # farad2; 

maxinh3 := d excit2; 

% Maxwell Lagrangian and energy-momentum current are assigned 
pform lmax4=4, maxenergy3(a)=3$ 
lmax4 := -(l/2)*f arad2~excit2; 

maxenergy3(-a) := e(-a) _ I lmax4 + (e(-a) _ I f arad2) ~excit2; 



% Be sure that you use a blank before the interior product sign! 
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